Short title: Baryons from u(3) 

Neutron to proton mass difference, parton distribution functions and bciryon 
resonances from dynamics on the Lie group u(3) 

01c L. Trinhamnier 
Department of Physics, Technical University of Denmark, 
Fysikvej, Building 307, 
DK-2800 Kongens Lyngby, Denmark, 
ale. trinhammer@fysik. dtu. dk 

We present a hamiltonian structure on the Lie group u(3) to describe the baryon spectrum. 
The ground state is identified with the proton. From this single fit wc calculate approximately 
the relative neutron to proton mass shift to within half a percentage of the experimental value. 
Prom the same fit we calculate the nucleon and delta resonance spectrum with correct grouping 
and no missing resonances. For specific spin eigenfunctions we calculate the delta to nucleon mass 
ratio to within one percent. Finally we derive parton distribution functions that compare well 
with those for the proton valence quarks. The distributions arc generated by projecting the proton 
state to space via the exterior derivative on u(3). We predict scarce neutral fiavour singlets which 
should be visible in neutron diffraction dissociation experiments or in invariant mass spectra of 
protons and negative pions in B-decays and in photoproduction on neutrons. The presence of such 
singlet states distinguishes experimentally the present model from the standard model as does the 
prediction of the neutron to proton mass splitting. Conceptually the Hamiltonian may describe an 
effective phenomenology or more radically describe interior dynamics implying quarks and gluons 
as projections from u(3) which we then call allospace. 



PACS numbers: 12.90.+b, 14.20.Gk, 14.20.Dh 

1 Introduction 

The quark flavour model has a missing resonance pro- 
blem since it predicts many more baryon TrA'^-resonanccs 
than observed [1, 2, 3]. The quark colour model (QCD) 
has a confinement problem to construct hadrons analy- 
tically from quarks and gluons. We are aware of the 
successes of QCD in pertubative domains [4, 5] and in 
lattice gauge theory [6, 7]. Nevertheless we want to stress 
the two first mentioned issues. The latter, the confine- 
ment problem, motivated a radical approach which solves 
the former as a by-product. We construct the dynam- 
ics in a compact space. The spectroscopy can live there 
and manifest itself in real space as different mass reso- 
nances. Then we will have confinement per construction 
and can hope to see quarks and gluons by projection 
from the compact configuration space to the laboratory 
space. We will present here a study of this idea. We shall 
call it the allospatial hypothesis from the greek word al- 
ios meaning another or different. In general terms the 
idea is that a hamiltonian description is the more nat- 
ural framework for spectroscopy and the lagrangian de- 
scription more suitable for scattering phenomena, see e.g. 
[8]. 

The defining equation is a group space Hamiltonian 
on u(3) where the toroidal degrees of freedom project 
out in colour quark fields and gluons come out by an 
adjoint projection. It turns out that there is also room 
for flavour in the model and we can reproduce an Okubo 
mass formula involving hypercharge and isospin. We note 
from the beginning that u(3) is not to be thought of as 



a symmetry group but as a conflguration space. 

The allospatial hypothesis is evolved in sect. 2. Here 
a speciflc Hamiltonian is stated, the underlying quanti- 
zations defined and a projection of the group space wave 
function into real space quark and gluon fields given. The 
projection is shown in appendix A to lead to fields trans- 
forming properly as respectively fundamental and adjoint 
representations of su(3). Appendix A concludes: "Prom 
the projection ... to laboratory space we recognize the 
toroidal generators as momentum operators. Thus when 
experimental production of resonances is of concern we 
see from space: The impact m,omentum generates the 
(abelian) maximal torus of the u(3) allospace. The mo- 
mentum operators act as introtangling generators. When 
decay, asymptotic freedom, fragmentation and confine- 
ment is of concern we see from allospace: The quark and 
gluon fields are projections of the vector fields induced, 
by the momentum form With this interpretation of 
the origin of quarks we see that the quarks are confined 
per construction since the Lie group of allospace is com- 
pact and thus cannot be projected in its global totallity. 
"You cannot peel an orange without breaking the skin", 
as differential topologists say, see fig. 1. 

The theory is unfolded in sect. 3. The Laplacian is 
parametrized in a polar decomposition analogous to the 
treatment of the hydrogen atom in polar coordinates. We 
shall see that the Laplacian contains a term naturally in- 
terpreted as a centrifugal potential. It includes the off- 
toroidal generators of the group and thereby carries the 
complexity needed to include spin and flavour spectral 
characteristics alongside the basic colour dynamics. As 
an example we calculate matrix elements of our Hamil- 
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tonian with symmetrized D-functions of specific angular 
momentum and find a very promising ratio when com- 
pared with TOA(1232)/'7lAr(939)- 

Energy eigenstates are discussed in sect. 4 and both 
approximate and exact solutions given for the neutral 
flavour N- and A-sectors together with tentative spin- 
parity assignments. The number and grouping of reso- 
nances agrees with all the certain (four star) resonances 
listed by the particle data group without a missing reso- 
nance problem. 

In sect. 5 we give a controversial interpretation rela- 
ting period doublings in the wavefunction to the creation 
of charge in the neutron decay. The mass shift related 
to this interpretation is quite promising with a relative 
mass difference of 0.13847 % between the neutron and 
proton mass predicted from the approximate solutions 
to be compared with the experimental value of 0.13784 
%. In sect. 6 we discuss experimental predictions. Apart 
from predicting all the observed, certain resonances our 
model predicts scarce neutral flavour singlets. In sect. 7 
we list some open questions together with selected parton 
distribution functions and in sect. 8 we give concluding 
remarks. 

The present work contains five appendicies. In ap- 
pendix A we give the projection to space. In appendix B 
we derive the spectrum of the off-toroidal generators of 
our group space Laplacian. In appendix C we describe 
the Rayleigh-Ritz method used to solve the exact case for 
neutral charge states and calculate analytically the ma- 
trix elements of the Rayleigh-Ritz method when this is 
used on trigonometric base functions. In appendix D we 
derive parton distribution functions. In appendix E we 
give a survey of those concepts in differential geometry 
most crucial for the calculations in appendicies A and D. 



2 The allospatial hypothesis 

The 'communication' between the interior dynamics 
and the spacetime dynamics runs through the exterior 
derivative of the state on the chosen configuration space. 
So the characteristics of the fields we wish to generate by 
the exterior derivative implies which configuration space 
should be hypothesized. We wish to generate projec- 
tion fields transforming under the SU(3) algebra with 
the fields possibly being electrically charged. This points 
to a configuration space containing both su(3) and u(l). 
Thus we choose the Lie group u(3) as configuration space 
and assume the following Hamiltonian 



he 
a 



(1) 



It is the hypothesis of the present work, that the eigen- 
states of (1) describe the baryon spectrum with u e u{3) 
being the configuration variable of a sole baryonic entity 
and a is a scale. Below we shall find exact solutions of 
(1) for alleged N-states and we shall discuss approximate 
solutions for both alleged N- and A-states. It should be 
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From algebra to group and vice versa 



FIG. 1: The algebra approximates the group in the neigh- 
bourhood of origo. The allospace in the present work is the 
Lie group u(3) in which the dynamics is supposed to describe 
the baryon spectrum as a set of stationary states of a certain 
Hamiltonian. Projecting from the group to the algebra gene- 
rates quark and gluon fields according to the representation 
space chosen. 



mentioned that, when unfolded, the structure of (1) car- 
ries degrees of freedom for both colour, spin, hypercharge 
and isospin. 

The Laplacian A in (1) is parametrized in the next 
section. The configuration variable u = e"^^, and the trace 
of its squared argument Trx^ — (P{e,u) is used in the 
representation independent potential. So the potential is 
half the squared geodetic distance [9] from the 'point' u 
to the 'origo' e 



d\e,u) = 9l + 9i + 9l 



-TT < 9. < TT, 



(2) 



where e'^^ are the eigenvalues of u. The geodetic distance 
may be seen as the euclidean measure folded into the 
group manifold [10]. The geodetic distance is invariant 
under translation in group space as it should be since the 
choice of origo is arbitrary 



d{v,uv) — d{e,v ^uv)=d{e,u). 



(3) 



Planck's constant h = 2'Kh enters the description via 
the canonical quantization behind eq. (1). The speed 
of light c and the scale a enter the description from a 
parametrization of the eigenangles 9 via either a time 
parameter t or a space parameter x — a9 = ct. This gives 
the scale huj = h9 = he/ a for the dimensionful eigenva- 
lues E. The ground state eigenvalue E = E/{he/a) of 
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FIG. 2; Periodic parametric potential (17) originating from 
the squared geodetic distance (2) from u to the origo e in the 
Lie group (Milnor [9]). The dashed curve corresponds to the 
Wilson analogue [11] of the Manton potential [12] and is not 
considered in the present work. 



It is via (6) and (8) that the physical dimensions enter 
the model (1). 

The state 5" projects out on laboratory space via 
d$ where $ — J^P and J is the Jacobian (13) of our 
parametrization. For instance the restriction of the mo- 
mentum form to the torus 

d$ = ^Pjd0j (9) 

generates a member of the fundamental representation of 
SU(3). The coefficients in (9) are the components of the 
alleged colour quark spacetime field whose transforma- 
tion properties are shown in appendix A. 

3 The theory unfolded 



the dimensionless edition (10) of (1) is of the order of 4.5 
predicting a promising scale a « 1 fm if the dimensionful 
ground state eigenvalue is identified with the nucleon E = 
939 MeV. Reversing the argument a is to be determined 
experimentally by fixing the ground state energy of (1) 
to that of the nucleon. It turns out that A = he/ a— 210 
MeV in fine agreement with the average QCD scale, e.g. 

4^5 -[13]. 

We denote the configuration variable u by a lower case 
letter and accordingly the Lie group space to make it 
clear that u is a configuration variable and not a priori 
a symmetry transformation. We parametrize the confi- 
guration variable by u = e*"'=-^'= where for A; = 4, . . . , 9 
are the off-toroidal Gell-Mann generators (14), whereas 
Tj for j — 1,2,3 are toroidal generators and aj — Oj. 

We have toroidal coordinate fields dj implied by the 
(abelian) toroidal generators 



d 



(4) 



The basic action-angle quantization condition {h = 1) 
with Tj — —id/dOj — —idjle reads 



(5) 



for z,j = 1,2,3, where dOi are the torus forms. 

The scale a enters the model when the three toroidal 
degrees of freedom are related to the three spatial degrees 
of freedom by the following projection 



Xj — aOj. 



(6) 



Thus we define parametric momentum operators pj pro- 
portional to the toroidal generators 



.,1 d h 



a do. 



J' 



(7) 



To study the stationary states of (1) let us rewrite the 
Schrodinger equation in a dimensionless form with E = 

E/K 



-^{u) = E^'(u). 



(10) 



The Laplacian from a polar decomposition on u(3) is [14] 



Ml 



Ssin^ i( 



(11) 



We recognize a term for toroidal kinetic energy and two 
terms which may be interpreted in the Hamiltonian as 
additional potentials: a constant global curvature [15] 
potential \ ■ (—2) = —1 and a centrifugal potential 



C 



1 ^ 



Kl 



Ml 



,2 1 



(0. 



(12) 



Further the van de Monde determinant, the 'Jacobian' of 
our parametrization is 



j=n2. 



1 



(13) 



and with h = 1 the off-diagonal Gell-Mann generators 
are [16, 17, 18, 19] 

K2 = X5, K3^X2, (14) 

Ml = Ae, M2 = A4, Ah = Ai. 

Kk commute as body fixed angular momentum opera- 
tors and Mfe 'connect' the algebra by commuting into 
the subspace of Kk 



whereby the canonical action-angle quantization of (5) 
translates into position-momentum quantization 



[aOi,Pj] 



{[a0i,a9j] = [pi,Pj] = 0) 



(8) 



[Mk,Mi] = [Kk,Ki] = -iK^ (cyclic in fc,/, to). (15) 

The presence of the components of K = {Ki, K2, K^) in 
the Laplacian operator opens for the description of spin. 
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Interpreting K as the allospatial spin operator is encour- 
aged by the body fixed signature of the commutation re- 
lations. The relation between space and allospace is like 
the relation in nuclear physics between fixed coordinate 
systems and intrinsic body fixed coordinate systems for 
the description of rotational degrees of freedom. Once 
and e.g. have been chosen to be independent and mu- 
tually commuting operators, both commuting with the 
Hamiltonian H of (1), the presence of M = {Mi,M2, A'h) 
in the Laplacian makes it possible to describe flavour via 
the spectrum of M^, see below. 

Substituting (11) into (10) we get the fully 
parametrized Schrodinger equation 



• 1 J^^J 



Kl + Ml 



+v{ei) + v{e2) + viOsJ'fiu) = E*(u), (16) 



where the potential ^cf spells out as a sum of periodic 
'chopped' harmonic oscillator potentials. The potential 
(17) is shown in fig. 2, 



v{9) = \i^- nnf^ S e [(2n - l)7r, (2n + l)7r], 
where n Cz Z. 



(17) 



We may assume the wavefunction to be an eigenstate 
of and and thus write it as a product of a toroidal 
part t(0) like the radial wavefunction for the hydrogen 
atom and an off-toroidal part T^m like the spherical 
harmonics Yim- With = {61,02,63) we thus write 

^'(m) = t{0i,92,9s) ■ TKMia4,a5,a6,a7,as,ag). (18) 

It should be noted that in determining the spectrum for 
wc cixploit the possibility of choosing specific eigen- 
states of hypercharge and isospin 3-component as shown 
in appendix B and thus instead of Tkm we might had 
written TkyIs ■ Multiplying (16) by the Jacobian J we 
introduce a new function 

$(w) = R{e) ■ Tkm where i?(6>) = J{d) ■ r(6>). (19) 

Then $ satisfies 



E 

i<j 



8sin=^ \{9i - 6j) 



= 2E$(m), 
where the 'euclidean' Laplacian is 



(20) 



(21) 



Now we integrate out the off diagonal degrees of free- 
dom (a4, as . . . , ag) to get an equation for the torodial 



part Rkm for specific eigenvalues K{K + 1) and of 

respectively and 

[-A, + v]RKM{.9i,e2,e3) = 2ERKM{6i,e2,e3). (22) 

Here the curvature and centrifugal terms of the group 
space Laplacian have been collected with the geodetic 
distance potential into a total potential V 



1 ^ 
V = -2 + -{K{K + l))+M^)Y, 



.^^.Ssin^\{6i-6i) 



+2{v{ei)+v{62)+v{ez)). 



(23) 



In the centrifugal term during the integration we 

have exploited the existence of the Haar measure over 
(a4, 05 ... , ag) together with the factorization (18). Fur- 
ther we have used that the off-toroidal part of the wave- 
function is an eigenstate of and together with the 
fact that the centrifugal term (12) is symmetric under in- 
terchange of the torus angles 9j. 

The centrifugal term leads to a mass formula of the 
well-known Okubo type [20] . In appendix B we derive the 
spectnun of + and show the following relation 
among quantum numbers 



K{K + 1)+M'^ = -[n + 



-3--y2-4ii. (24) 



Here y is the hypercharge, is the three-component of 
isospin and n is an integer which we may call hyperdi- 

mension. It is natural in the present framework to classify 
the eigenstates according to the three independent values 
of n, y and i^. However we can make a transformation of 
this classification into the familiar one by rewriting the 
expression (24) and choose the sum of hyperdimension 
and hypercharge to be a constant. For n + y = 2, which 
yields the lowest possible K{K -h 1) + M^, we get 



i^(ir+l)+M2 = |+(fc|+m^)-|y+4 



j(i + l) 



(25) 

Since {K^ + M|) commutes with both Y and we get 
for a given value of (A:| -|- m\) 



K{K + 1)+M'^ = a' + b'y + c' 



y2 - i{i + 1) 



. (26) 



Equation (26) is the famous Okubo mass formula that 
reproduces the Gell-Mann, Okubo, Ne'eman mass rela- 
tions within the baryon N-octet and A-decuplet [17, 20, 
21, 22] independently of the values of a', b' and c'. Of 
course this is only so if one chooses the same toroidal 
wavefunction for all members of a given multiplet. In 
practice the SU(3) symmetry breaking in (26) will be 
influcncc;d by the 0-dcpendcncc in (23) because different 
values of K^-|- lead to different values of the centrifu- 
gal potential and thereby influence which span of toroidal 
energy eigenstates will project out on a specific angular 
momentum eigenstate in the laboratory. 
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3.2 Total angular momentum 

All members of an SU(3) multiplet have the same total 
angular momentum and therefore it would seem natural 
to let this be carried by the structure of the toroidal wave- 
function. As a possible base one may try ZJ-functions. 
The /^-functions are eigenfunctions of angular momen- 
tum [23] 



(27) 



Now the D-functions are functions of the Euler angles 
describing rotations between a fixed coordinate system 
and body fixed coordinates [24]. They are simultane- 
ously eigenfunctions of the three-component of angular 
momentum in coordinate space and in body fixed space 
with eigenvalues —m and —k respectively, where m and 
k may vary in integer steps from —j to j [25, 26]. De- 
spite the notation no a priori relation is intended between 
k, m and K, M. In a specific example below however, 
we shall test the identification j = K. In the present 
connection we may interpret the fixed coordinate system 
as the laboratory space and the body fixed space as al- 
lospace. Since a priori there is no prevailing direction in 
laboratory space the states should have no specific three- 
component of angular momentum. Thus states should be 
constructed by summing over m. And since the labeling 
of the eigenangles in allospace is arbitrary, we should ex- 
pand on states that are symmetric under interchange of 
the eigenangles. Therefore a possible expression for the 
nucleon ground state wave function with j = k = 1/2 
could be 



^1/2^(^1' ^2, ^3) — (^1/2(^1' ^2 



1) 



+tl$2(^3,ei,e2))M/2, 



(28) 
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FIG. 3; AH observed certain neutral flavour baryons (boxes) 
compared with approximate allospatial predictions (lines) 
from eq. (31) with label combinations based on table 2. The 
dashed line is a predicted neutral charge state without charged 
partners. The red lines mark states with augmented contribu- 
tion in level 3. The boxes indicate the experimental range of 
peak values [27], not the widths which are much larger. The 
systematic errors in the approximation are of the order of 50 
MeV. Digits at selected allospatial states are the toroidal la- 
bels m, n. Note the fine agreement in the grouping and the 
number of resonances in both sectors with just one fitting at 
p= l',2',3. 



,1/2 



IS a sum 



where Ni^2 is a normalization factor and t^^^ 
over the two Z?- functions for the two possible values 1 /2 
and -1/2 of m corresponding to total angular momentum 
J = 1/2 



1/2 
1/2 



D 



1/2 

1/2,1/2 



D 



1/2 

1/2,-1/2 



oce-'«^/2cos(^)e-'^3/2 
^ 2 ^ 



_e-^«i/2sin(^)e*«-^/2. 
^ 2 ' 



(29) 



< JT, 



1/2 
1/2 



Hr I JT^j^ > in (30) is 5.93 and fitting 



this to the nucleon rest energy 939.6 MeV corresponds 
to a scale a « 1.2 fm. It is also reassuring to note 
that the numerical eigenvalue 5.93 is rather close to the 
value 0.5111/0.0855 = 5.98 of the nucleon ground state 
in a favoured choice of parameters for lattice calculations 
shown in table 5 of [7]. 



For total angular momentum j = 3/2 there will be 
four terms in (29) corresponding to the allowed values 



of m 



3 1 



It is encouraging for this line of 



investigation to note that 



< JT. 



3/2 



3/2 



< JT, 



1/2 
1/2 



Hr 
Hr 



"'^3/2 
•^^1/2 



= 1.32 



™A(1232) 
"^Ar(939) 



= 1.31. 



(30) 

Here H^; is the allospatial Hamiltonian in the exact form 
in (22) and the result in (30) is for the lowest possi- 
ble value of K^-l- hke in table 1. The value of 



3.3 Energy eigenstates 

Equation (22) can be solved by a Rayleigh-Ritz 
method, for details see appendix C. The eigenvalues of 
the full equation (22) are listed in table 1 for the low- 
est possible value of K^-|- M^. Here we have expanded 
R{0) on Slater determinants of trigonometric functions 
of increasing order (44). In table 4 consistency is shown 
by comparison with exact results for (22) by expansions 
on Slater determinants of parametric eigenfunctions (33). 
In table 1 are shown also possible identifications with 



baryon resonances. We see that all the N-resonances 
which the Particle Data Group considers certain [27] are 
present except N(1520). Thus the general tendency is 
promising but the predicted energy values are in need 
of improvement. Such an improvement is seen if instead 
of comparing the energy of the electrically neutral part- 
ners of each N-resonance one looks at the energy of the 
charged partner. The energies of the charged partners 
are calculated below for an approximate case and the 
comparison is shown in fig. 3. In sect. 6 we shall return 
to the question of N(1520). 

4 Approximate solutions 

Here we omit the constant curvature potential and the 
centrifugal potential. They are of the same order of mag- 
nitude and of opposite sign, see sect. 5. We may then 
expect to keep the general structure of the spectrum. We 
shall introduce a new degree of freedom to be related to 
the electric charge. 

Having integrated out the off-toroidal degrees of free- 
dom in (19), we are lead to an approximate edition of 
(22), namely the following seperable equation 




2 ^ de'i 



R{9i,92, 03) 



ER{9i 



3). 



(31) 



With both the Laplacian and the potential now being 
just sums, the equation separates into three Schrodinger 
equations for parametric eigenfunctions with periodic 
boundary conditions 



(t)i{9j) = ei(j)i{9j), 



-TT < 9, < TT. 



(32) 

Figure 4 shows to the left the first eight parametric 
eigenfunctions. The eigenstate ^' in (18) is symme- 
tric in the eigenangles as the labelling of angles in the 
parametrization is arbitrary. J is antisymmetric. In to- 
tal R{9i, 92,93) is antisymmetric and may be constructed 
[28] as a Slater determinant of three freely chosen ortho- 
gonal {I < m < n) parametric eigenfunctions. 



h) 



h) M03) 

92) 0m(^3) 



M^l) 4>n{0i) 



(33) 



Substituting the iiTM-indices each state thus gets three 
toroidal labels I, m, n, the labels of the parametric eigen- 
functions. The eigenvalues E of (31) are sums of three 
different parametric eigenvalues E — J2i=imn^i- 
instance we will have E — ei + 62 + for 



Tl23{0) = - 



h{9i) 01 (03) 

h{9i) Md2) M^s) 

h{9l) MO2) 03(^3) 



(34) 



FIG. 4: Periodic parametric eigenfunctions (32) from which 
the toroidal part (34) of the wavefunction of approximate 
states are constructed. Eigenstates of diminished eigenval- 
ues are shown to the right. The exact functions in (22) can 
be found from an expansion on such states or on a base set 
(44) mimicking these parametric eigenfunctions by sines and 
cosines of increasing order. 



We may interpret the antisymmetry under interchange 
of the columns in the Slater determinant i? as a colour 
antisymmetry under interchange of the eigenangles of u. 
To keep a uniform probability measure in group space 
when seen from parameter space, it is $ that we would 
like to project to space and not Because the external 
derivative is linear d$ inherits the antisymmetry of 
and carries it through to the projection space, e.g. 



4.2 Parity 



(35) 



We might expect a 'winding number' effect in txm due 
to the increasing number of oscillations in the base func- 
tions. Such a winding number effect may be expressed 
in specific partial waves to which the allospatial state re- 
lates in space. This would explain the general trend of 
increasing total angular momentum with increasing en- 
ergy. The idea is supported by table 5 where we see the 
allospatial states together with spin and parity of the 
N-resonance candidates. We can give an ad hoc assign- 
ment of parity according to each particular combination 



of toroidal labels l^m,n, namely 

P= (-1)™-' • (-1)"-™, (36) 

which in general reduces to P = (—1)"^'. The assign- 
ment (36) gives a rather consistent correlation to the 
observed resonances with only a few exceptions namely 
for the singlet states like 135 (see below) which seems 
to mix with ordinary doublet N-states to yield a state 
with reversed parity. There also seems to be a correla- 
tion between increments in toroidal label number and in- 
crements in angular momentum, a tentative trend being 
A J = Am+An in series like 125, 127, 147 |-H 
and 134,136,156 S". I"- This confirms the 
winding number effect. Using the parity assignment (36) 
and the trend AJ = Am we predict from the series 
235,237,(257) ^ i-,|-,(^-) a value of ^- for 
the allospatial state 257 with approximate energy 2542 
MeV in accordance with the three star rated N(2600) 
[27]. Although promising the spin-parity assignments in 
table 5 is just a tentative number game. We shall now 
give a more specific suggestion. 

We return to the projection a9j = Xj and construct 
parity eigenstates 

T+{e)=T{e) + T{-e) a.ndT-{e)=T{e)-T{-e) (37) 

based on states like (28). Identifying the operation of the 
parity operator P on x with a similar operation on the 
allospace eigenangles 6 would ensure that the states 
have intrinsic parity ±1 

P(r±(6i)) = T±(-e') = ±T±(6i). (38) 

Since [P, H] = states are possible with simultaneous 

1 /2 

P and E eigenvalues. A state like T^j2 (2^) then 
split into two with specific parity. The presence of two 
different parity states for each angular momentum fits 
very well with the pattern of the observed certain N- 
resonances. Such states might be the spin-parity-definite 
states through which the spin-parity-indeterminate en- 
ergy eigenstates of parametric function Slater determi- 
nants like (33) communicate with the partial wave am- 
plitudes of the experimental analysis. 



4.3 Bloch states 

Now let us return to the question of solving (31) and 
(32). The basic trick in exploring the full structure of 
m(3) and solving the parametric eq. (32) is to introduce 
Bloch waves (/)(0) = e^'^^g{9) and get a set of eigenvalue 
bands, see fig. 5. Here g is the 2tt periodic part belonging 
to the Bloch wave number k. Choosing k = gives us 2tt 
periodicity of the parametric functions corresponding to 
the periodicity of the parametric potential. For the state 
with lowest eigenvalue in this case E — ei + 62 + £3 — 
4.474172455, see table 6. If we choose K2' = ±i in level 




FIG. 5: Reduced zone scheme [29] for allospatial parametric 
eigenvalues. The black dots represent the values for the un- 
stable neutron state (left) and the proton state (right). For 
clarity the variation of the eigenvalues with wave number is 
grossly exaggerated for the lowest states. 



2 it will lead to a lower lying parametric eigenstate which 
we can couple to a 'compensating' Ky = in level 1. 
We thus couple two Bloch waves (f)v{0) = e^'^^^gi'{9) and 
4'2'{(^) — e^'^^^ 92' {(^) from the two lowest lying eigenvalue 
bands in the reduced zone scheme to the right in fig. 5. 
This combination is possible because the energy gained 
by the period doubling in level 2 is larger than the energy 
lift needed to cause period doubling in level 1. Note 
that the parametric eigenstates of the first and second 
levels now have 47r periodicity in the eigenangles. This 
is possible since we only require ^'■^ and not nescessarily 
\I' to be single valued on u(3) after a single "rotation" in 
the parameters 9j (cf. the case of a spin ^ particle like 
the free electron which needs a Air rotation for its state to 
return to its original value [30]). Other values of k than 
and ±^ are not in accordance with the single valuedness 
of 'I'^. The lower eigenvalue after this coupled decay is 
E" = 61+62+63 = 4.467985519 which we interpret as the 
proton ground state. With this interpretation we predict 
from our approximate calculation the relative difference 
between the neutron mass m„ and the proton mass rup 

m^-mp _ E- E" _ 4.474172455 - 4.467985519 
^p ~ E" ^ 4.467985519 

== 0.138473%. (39) 

The result in (39) is in close agreement with the experi- 
mentally determined value [27] 

run-mp _ 939.565346 MeV - 938.272013 MeV 
TTip ^ 938.272013 MeV 

= 0.137842%. (40) 
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It is encouraging that the approximate calculation gives 
a too high mass splitting. One may namely expect the 
eigenvalue bands to narrow slightly in an exact solution 
because of the spikes from the centrifugal potential which 
has been omitted in our above calculation and which will 
effectively leave the levels in a deeper potential well. On 
the other hand one may fear that the resulting level shifts 
might be too large in order to keep the fine agreement 
indicated by (39) and (40). However, the neutron state 
shifts only two per cent downwards to 4.38 for the ex- 
act solution in table 1. Actually the Jacobian (13) 'kills' 
the wave function exactly at the location of the singu- 
larities of the centrifugal potential, so we may hope that 
the shifting effect on the proton state will be just what 
we would need in order to diminish the still significant 
discrepancy between the prediction in (39) and the very 
accurate observations in (40) where the uncertainties on 
both masses are ±0.000023 MeV which means that the 
uncertainty lies on the last digit in (40). 



5 Charge interpretation and systematic error 

We interpret the period doublings in the diminished 
ground state topology as 'creator' of the proton charge 
in the neutron decay. Thus k, is interpreted as a charge 
quantum number degree of freedom which express the 
coupling between allospatial structure and the spacetime 
fields of the electroweak interaction. The interpretation 
is supported by the structure of the implied proton state 



''l'2'3 



e-'"-''gv{ei) e-'--'-gv{92) e-'^'^grie^) 



l) 



l) 



2) 



3) 



3) 



(41) 



We see that the period doublings induce u(l) factors 
curled up in the Slater determinant thereby "trapping" 
the charge topologically. Note that u(l) is the gauge 
group of the electromagnetic interaction. The diminished 
parametric function in the second row of the determinant 
in (41) is shown to the right in fig. 4 and the lifted one 
in the first row is shown in fig. 6. 

The decay seems to couple inherently to spin degrees 
of freedom since the two Bloch wave phase factors render 
the measure scaled wave function 



l) e^"^'^g2'{e2) 6^5^332,(03) 



, (42) 



with a structure resembling sums of 13- functions [32] . For 
instance 

d\ ,{91,62,9^) = 

2 ■■ 2 

e-'^'^'d\ _ 1 {92)e'^^'' = -e-*59i sin(^02)e'5''^ (43) 
2 ' 2 2 




FIG. 6: The first lifted parametric eigenstate with periodicity 
47r. The eigenvalue is found by collocation with 1500 points 
in Matlab and the corresponding eigenstate is generated by 
solving in Maple the one-dimensional Schrodinger equation 
(32) for the particular eigenvalue. See [31] for programs. 



is related to a particular rotation of spinors for the case 
of spin i [33, 34]. We may even imagine that the decay 
with its simultaneous period doublings in two parametric 
levels generates two component structures to carry away 
the excess energy. It is beyond the scope of the present 
work to show that these structures might be the electron 
and its antineutrino [35]. 

Instead of expansions on the parametric eigenfunctions 
(33) we could also use the Rayleigh-Ritz method for solv- 
ing the full (22) and the approximate (31) eqs. by expan- 
sions on functions of the type 



Upqr\9i, 6*2, 



cosp^i COSp6'2 cosp03 
sin q9i sin q92 sin qO^ 
cos r9i cos r02 cos r9^ 



lijk cos p9i sin q9j cos r9k , 



(44) 



where p 0, 1, 2, . . . ; g = 1, 2, . . . ; r = p + l,p + 2, . . . ; 
and similar with 0, 2 or 3 sines. The details are laid out 
in appendix C and it turns out that all matrix elements 
for this set of base functions can be calculated analyti- 
cally. The functions (44) have a structure quite similar 
to the Slater determinant (33). The increasing number 
of oscillations in the parametric eigenfunctions of (33) is 
modeled by the increasing order of the cosines and sines. 

The order of magnitude of the omitted terms in (31), 
i.e. the constant curvature potential and the centrifugal 
potential, may be evaluated from these base functions. 
For instance 



< boil I C 1 6011 > / < boil 1 6011 > ^ 



(45) 



where C is the centrifugal potential for — ^(^-1-1) and 
JVP = ^ like in table 1. The centrifugal potential more 
or less cancels the curvature potential -1 and the result 
-1/3 of the two terms taken together is small compared to 
the geodetic potential which yields 7r^/2, see (C12) and 
(C9). 



9 



6 Experimental predictions 

We determine the scale A for the approximate solu- 
tions from the proton rest energy K = E/Yj = 938.3 
MeV/4.468 = 210 MeV. Based on this scale the pre- 
dictions are shown in fig. 3. The allospatial spectrum 
agrees with the number and grouping of all the certain 
resonances in the N- and A-sector. It should be stressed 
that all the observed certain resonances are accounted 
for by allospatial states. By 'certain' we mean all the 
well established resonances with four stars in the particle 
data group listings [27]. Relative to the exact solutions 
in table 1 the omissions in the approximate solutions re- 
present systematic errors of the order of 50 MeV. 

The relative mass splitting between the neutron and 
the proton resulting from the period doublings in (39) is 
0.13847 % to be compared with the experimental value 
of 0.13784 %. The slightly lighter (as opposed to heav- 
ier) mass of the charged proton is natural in the present 
hypothesis. States with two even labels like 1,2,4 have 
the possibility of two parametric eigenvalue diminishing 
period doublings. This will give states with a multiplicity 
of four. The state 1, 2, 4 is interpreted as the A-resonancc 
with two diminutions for A++. Again the charged part- 
ners are naturally the lighter ones. The state 2, 4, 6 is 
the first example where all three parametric states can 
undergo diminution induced period doublings. But the 
complex exponentials are now common factors in each 
column of the Slater determinant so they factorize out 
thereby releasing the topological trap. The determinant 
regains a parametric 2it periodicity and the state may 
eventually contribute to neutral resonances. 

The grouping of the resonances agrees quite well 
with the observed certain states. Only one certain N- 
resonance is missing in the group of three resonances in 
the domain around 1500 MeV. However the approximate 
treatment in (31) predicts a neutral singlet 1,3,5 at 
1510 MeV exactly in that area. The exact treatment in 
(22) places the singlet at 1526 MeV, see table 3. This 
state is thought to mix with the other two N-resonances 
to give the total of three N-resonances in the group. The 
next singlet state 1,3,7 is predicted in the 'desert' area 
between 1700 MeV and 2100 MeV. In the approximate 
case of (31) the resonance comes out at 1965 MeV and 
in the exact case of (22) it comes out at 2051 MeV 
(table 3). No certain N-resonance is observed in this 
domain. Being close to the observed resonance domain 



around 2100-2200 MeV the state 1,3,7 may hide itself 
by mixing with ordinary N-states. Thus it could explain 
the new resonance N(2040) in m^^- invariant mass 
spectra from J/* — >■ p'K~n [36]. The singlet 5,7,9 lies 
just above the free charm threshold for baryonic decay 
into S+(2455)D- at 4324 MeV (2454+1870 [27]) and 
should be visible (together with 3,5,11 at 4652 MeV) in 
neutron diffraction dissociation experiments like those 
in ref. [37]. It should also be visible in 'yn — > p7r~ 
photoproduction experiments like [38]. Other, lower 
lying singlets shown in table 3 might be visible in 
invariant mass m^^- from B decay experiments like [39, 
40]. 



7 Open questions 

The present work has three open ends: a very specific, 
momentarily halted end; a far horizon perspective and 
finally a more immediate and accessible opening. 

1. It has not been possible to solve the general eq. (22) 
for the alleged electrically charged states. A suitable set 
of base functions analogues to the uncharged case (44) 
has not been found. 

2. The introduction of the charge degree of freedom in 
(41) is supported by the fact that u{l) is the gauge group 
generating the electromagnetic interaction. This point is 
in need of elaboration. 

3. The projection mediated by the exterior derivative 
relates toroidal parameters to spacetime parameters and 
transfers structure from the group configuration space 
into transformation properties of the spacetime fields. 
However we also need to understand how the topology 
of the allospatial state transfers into angular distribu- 
tions in spacetime partial waves. We expect _D-functions 
to be involved. They are representations of the rota- 
tion group o(3) which happens to be doubly covered by 
su(2) of which there are three 'copies' in su(3), e.g. re- 
presented by the three spins U , V and /. We may then 
assume that the spin structure in allospace is excited via 
the generators of the rotation group in laboratory space. 
If we apply the so-called Clebsh-Gordan series [41] twice 
we can reduce a triple product of _D-functions to a sum 
over f-functions presumably pertaining to the toroidal 
wavefunctions with specific angular momentum j 



j' 

Y^C{hj2j'\ kik2)C{jij2j'; mima) Y^Cifjaj; fci-hfca, kz)C{j'jzj\ rm+ma, m3)I>^^+fe^+fe3_„^+„^+„3. (46) 
j' i 
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A systematic investigation of total spin and parity 
eigenstates based on Z?-functions like (37) seems an ob- 
vious task to undertake from here thus attacking point 
3, but we would like to sort out the relation between the 
structure of the D-functions and the toroidal labeling of 
states, on which our present interpretation of charge mul- 
tiplets rests. After all the energy eigenstates constructed 
from expansions on parametric function Slater determi- 
nants agree surprisingly well with the observed states as 
seen in fig. 3. Further the toroidal parton distributions 
traced out in appendix D and shown in fig. 7 compares 
rather well in shape with the proton valence quark dis- 
tributions in figs. 16.4 of [42] and 6.14 of [43]. 




FIG. 7: Parton distribution functions for a first order ap- 
proximation to an alleged proton ground state of (1). The 
distributions are traced out by specific toroidal generators 
Tu = (|, 0, — 1) and Td = (— |, 0, — 1) via the exterior deriva- 
tive (9), see appendix D. The curves generated by Tu (green 
solid) and Td (red dashed) compare both in peak value and 
relative integral content with the proton valence quark distri- 
butions Uv and d„ shown in the insert from [42]. 



7.2 On Interaction 

Ansatz 

Baryons interact via the harmonic potential, their 'dis- 
tance' being measured by the length d{u,u') of the 'rela- 
tive' geodetic where u and u' are the configuration vari- 
ables of the respective baryons. This leads to two con- 
jectures. 

Conjecture 1 

The deuteron should be the ground state of the follow- 
ing Schrodinger equation 



a 



(47) 



The isotopic landscape should in principle be derivable 
from a Hamiltonian something like 



He 
a 



Ha — Hz+N — 

Z,N 

^2 



A{A-l)/2 



k<l 



(48) 



Conjecture 1 should be tractable, at least numerically, 
by exploiting the translational invariance of the distance 
measure. Here one would rewrite 



d{u,u') — d{e,u ^u') — d{e,u^u'), 



(49) 



where in the last expression we have used the unitarity 
of M. In that way we end up with a system of six differ- 
ential equations in the two sets of eigenangles. This set 
is coupled by the potential term 



1 



d2(e,u^u') 



2 



(50) 



Conjeeture 2 



Conjecture 2 may be tested by studying the decay of tri- 
tium to helium-3. According to our allospatial hypothesis 
this decay involves period doublings analogous to the de- 
cay of the neutron and thus in principle one should be 
able to calculate the relative mass difference of f H and 
fHe without knowing the actual mechanism behind the 
creation of the electron and its antineutrino. 



8 Conclusions 

We have presented a radical approach to the confine- 
ment problem by shifting the dynamics of strong inter- 
action baryon spectroscopy to a compact configuration 
space. We have developed a hamiltonian framework in 
which to formulate this dynamics and suggested how the 
dynamics communicate with spacetime kinematics. We 
have called it the allospatial hypothesis. The model has 
no fitting parameters, except the scale A which turns out 
to be 210 MeV corresponding to a length scale a 1 fm. 
The number of resonances and their grouping in the N 
and A sectors are reproduced well by eigenstates of a 
Schrodinger equation on u(3), called allospace. A quite 
accurate prediction of the relative neutron to proton mass 
shift 0.138 % follows from approximate solutions to the 
same Schrodinger equation. A projection of states to 
space is given via the exterior derivative. This projec- 
tion has shown to yield parton distribution functions 
that compare rather well with those of the proton va- 
lence quark distributions already in a first order approx- 
imation. A kinematic parametrization for the projection 
gives a natural transition between a confinement domain 
where the dynamics unfolds in the global group space and 
an asymptotic free domain where the algebra approxi- 
mates the group. A promising ratio between the A(1232) 
and N(939) masses has been calculated based on specific 
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D-functions. We expect the allospatial energy eigenstate 
spectrum to project into partial wave amplitude reso- 
nances of specific spin and parity via expansions on spe- 
cific combinations of D-functions. Singlet neutral flavour 
resonances are predicted above the free charm thresh- 
old of S+(2455)Z)~. The presence of such singlet states 
distinguishes experimentally the present model from the 
standard quark model as does the prediction of the neu- 
tron to proton mass splitting. 

Acknowledgments 

We thank Torben Amtrup, Jeppe Dyre, Jakob Bohr, 

Mogens Stibius Jensen, Mads Hammerich, Povl Holm, 
Hans Bruun Nielsen, Holger Bech Nielsen, Norbert 
Kaiser, Vladimir B. Kopeliovich, Sven Bj0rnholm and 
Bo-Sture Skagerstam for key comments. And we thank 
Erik Both for the meticulous transformation of our 
manuscript into I^I^t;X. We have also had help from Geof- 
frey C. Oades, Abel Miranda, Poul Werner Nielsen, Jens 
Hugger, Bo Gervang, Ivan Stich, Karen Ter-Martirosyan, 
Yurii Makeenko, Dmitri Boulatov, Jaime Vilatc, Pedro 
Bucido, Manfried Faber, A. Di Giacomo, Gestur Olafs- 
son, M. Blazek, P. Filip, S. Olejnik, M. Nagy, A. Nogova, 
Peter Presnajder, Vlado Ccrny, Juraj Bohacik, Roman 
Lietava, Miroslava Smrcinova, Hans Plesner Jacobsen, 
Vracheslav P. Spiridonov, Norbert Kaiser, Poul Olesen, 
Ben Mottelson, Andreas Wirzba, Niels Kja^r Nielsen, Bo- 
Sture Skagerstam, Per Salomonson, Palo Vaslo, Lissi Re- 
gin, Berit Bj0rnow, Elsebeth Obbekjaer Petersen, Hans 
Madsb0ll, Karl Moesgen, Victor F. Weisskopf, J0rgen 
Kalckar, Leo Bresson, Knud Fjeldsted, Jens Bak, Mo- 
gens Hansen, Jane Hvolba^k Nielsen, Jens Lyng Petersen, 
Peter Hansen, J0rn Dines Hansen. 

Appendix A: Projection to space 

For each element u G u(3) we have a corresponding left 
translation lu on v € u{3) 

lu{v) = uv, (Al) 

and for any left-invariant vector field X we have [44] 

= d(Z„)„(X„). (A2) 

In particular we have for our toroidal coordinate fields 
when comparing with (4) 

dj \u-e= d{lu)e{dj |e) = udj |e . (A3) 

Thus the exterior derivative d acts as the identity on left 
translations at the origo e. We now expand the exterior 
derivative of the measure scaled wave function $ = Jip 
on the torus forms (5), i.e. 

d^ = i)jdej, (A4) 



where the coefficients are the local partial derivatives [45] 
^l;i{u) = d^u{dj), i = 1,2,3. (A5) 
For the coefficients we have by left invariance (A3) 

^j{u) = d^uidj) = [$] = ud,u m 

= ud^eidj) = uipj{e). (A6) 

The sum of the torus form partial derivative coefficients 
will inherit the left invariance 

V'(m) = V'l(w) -I- 1p2{u) + Ipsiu) 

= u(V'i(e) -I- V2(e) -I- V3(e)) = wpie). (A7) 

Now in particular ipj{e) = d^e{dj) = d^/dOj 
belongs to the tangent space TMg of the max- 
imal torus M at e and therefore so does 
their sum ip{e) as in general i^jiu) G TM^. 
The set of generators {iTj} are the coordinate field 
generators dj which also constitute an induced base 
from parameter space 

d 

dj\u= Qf\n= d{exp)^^p-^^){cj), (A8) 

where {c, } is a set of base vectors for the parameter space 
for the torus, see fig. 9, appendix E. In our interpretation 
we identify {cj} with a base for a fundamental represen- 
tation space for the colour algebra SU(3). We may thus 
introduce complex valued components ipj for the colour 
vector tjj and write 

_ . _ m 

tp = V'lCi -I- V2C2 + V3C3 = < V'2 > • (A9) 

Us J 

In the above representation u will be represented by a 
3x3 matrix U. For rotations under V G SU{3) we then 
have 

c, ^ dj = Vcj (AlO) 

and 

U ^U' = VUV-\ (All) 

From (A7), (AlO) and (All) we can derive the transfor- 
mation property of tp{u). 

VUV-'^Vjj{e) = VUjj{e) = F^(u), (A12) 

which shows that the differential coefficient vector tp 
transforms as a colour vector in the fundamental re- 
presentation under SU(3) rotations. In other words left 
invariance in group space projects out as SU(3) rotation 
in projection space. We thus interpret ijj in (A7) as a 
quark field with three colour components which may be 
projected on a specific base like in (A9). 
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Likewise the gluon fields may be seen as resulting from 
a projection on adjoint representation spaces of an expan- 
sion of the momentum form corresponding to the full set 
of eight generators Xk needed to parametrize the general 
group element u = 6*"'=^'= - separating out a phase factor 
[46]. Thus for each generator we have left invariant 
vector fields da defined as 

d 



da = 



,=0= uiTa, 



(A13) 



where = —id/daa = —ida\e- We now choose the set 
{Ta} as a base for the adjoint representation. This base 
transforms under V G SU{3) hke 



T' = VTaV- 



(A14) 



Analogous to (A4) we expand the external derivative d$ 
on forms related to the left invariant vector fields da to 
get the adjoint projection field 



A{u) = Y,d^u{da). 



(A15) 



We want to show that A transforms according to the 
adjoint representation. First we have the equivalent of 
(A7) 

A{U) =J2d^u{da) =Yl9a\u m 

a a 

= 53 Udale [*] = U ^ d^eida) = uA{e). (A16) 
a a 

Here we understand in analogy with (A9) that 

A{e) = Aa{e)Ta, (A17) 

where again Aa are complex valued components. We may 
then proceed to show the adjoint transformation property 

of 

A{u) A{u'y = U'A{e)' = [/'ia(e)T^ 
= VUV-^Aa{e)VTaV-^ = VUV-WAaie)TaV-^ 

= VUA{e)V-'^ = VA{u)V-'^, (A18) 

which corresponds to the gauge group rotation transfor- 
mation property of the gluon fields B [47] 



b; = vB^v-' 



:id,v)v-\ 



where 



V 



-iaa{x)Ta, 



(A19) 



(A20) 



We note that as space time fields the gauge fields also ac- 
quire a term representing the variation along spacctime 
translations as represented by the second term in (A19). 
Note further that translational invariance in group space 
corresponds to an SU(3) rotational invariance in repre- 
sentation space and thereby the translational invariance 
of the interaction potential (3) in group space through 
the projections (6) and (9) reflects the gauge invariance 
of the fields in laboratory space. 



1. Kinematics 

Using (6) and (7) and suppressing the centrifugal and 
curvature terms like in (31), the projected Schrodinger 
equation reads 



1 2 9 d 

2 dxj dxj 



Ea 



he 



R{yi/a) = — E(x/a). 

(A21) 

In the high energy limit, taking S ^ oo to be the la- 
boratory energy and letting a ^ while the allospatial 
eigenvalue E = ^ remains fixed, we see that the pro- 
jected potential blows up and the kinetic term cools down 
to restrict the support of R to the neighbourhood of origo 
where the algebra approximates the group, see fig. 1. 
Compactification loosens up. This is a sign of asymptotic 
freedom. In the other extreme taking -> i?o = TOqC^ 
we will have a increasing, thereby lowering the poten- 
tial such that the wavefunction spreads out in all of the 
group. Adjusting E in the experimental set-up auto- 
matically determines whether one sees the deglobalized 
(asymptotically free) or the globalized (confining) char- 
acteristics prevail. 

Prom the projection (6) and (7) to laboratory space 
we recognize the toroidal generators as momentum oper- 
ators. Thus when experimental production of resonances 
is of concern we see from space: The impaet momen- 
tum generates the (abelian) maximal torus of the u(3) 
allospace. The momentum operators act as introtangling 
generators. When decay, asymptotic freedom, fragmen- 
tation and confinement is of concern we see from al- 
lospace: The quark and gluon fields are projections of 
the vector fields induced by the momentum form 



Appendix B: The spectrum of and 

In the allospatial interpretation the presence 
of the components of K = (i^i, i^a) and 
M — {Mi,M2,M^) in the Laplacian opens for the 
inclusion of spin and non-neutral fiavour. It can be 
shown [14] that the components commute with the 
Laplacian as they should since the Laplacian is a 
Casimir operator and therefore the generators of the 
group commutes with it. They also commute with the 
geodetic potential, so 



[iffc,H] = [Mfe,H] =0, 



(Bl) 



where H is the Hamiltonian in (1). Further 

[Kk,\^^] = [KkM^] = [Mfc,M2] = [Mfc,K2] = 0. (B2) 

Thus we may choose K^, /Cs, as a set of mutually 
commuting generators which commute with the Hamilto- 
nian H. To solve the general problem (1) or (10) we need 
to specify the spectrum of the involved operators. The 
well-known eigenvalues of and Kj, follow [48] from the 
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commutation relations (15). Here we choose to interpret 
K as an interior angular momentum operator and allow 
for half- integer eigenvalues. Instead of choosing eigenval- 
ues of K3 we may choose I3, the isospin 3-component. 

To determine the spectrum for we introduce a 
canonical body fixed 'coordinate' representation [49] 



Ki = a62P3 — a03P2 = hXr 
K2 = a6iP3 - aO^pi = hXs 
K3 = a6\P2 — a92Pi = h\2. 



(B3) 



The remaining Gell-Mann generators are traditionally 
collected into a quadrupole moment tensor Q, but we 
need to distinguish between the two diagonal components 

^ {el + el-20l) + -^'^{pl+pl-2pl) = M, 



2^v"i -2 2^3 

Q^lh=\{el-6l) + Y^^{pl-pl) = X3 



(B4) 



and the three off-diagonal components which we have 
collected into M 



Ms/n = 6*1612 + pPiP2 = Ai 
a? 

M2/n = 6361 + j^P3Pi = A4 



Mi/n = 6263+ j^P2P3 = Ae- 



(B5) 



M is a kind of Runge-Lcnz 'vector' [50] of our problem. 
This is felt already in its commutation relations (15). 
We shall see in the end (B13) that conservation of 
corresponds to conservation of particular combinations 
of hypercharge and isospin. 

For the spectrum in projection space we calculate the 
SU(3) Casimir operator using the commutation relations 
(8) 

C, = ^(K' + M' + Ql + Ql) = -3 + 3^(2He)2, (B6) 

where the Hamiltonian He of the euclidean harmonic os- 
cillator is given by 



2He = — p 

n 



2 I 



with p = (pi,P2,P3), d = {61,62,93 



(B7) 



We define canonical annihilation and creation opera- 
tors 

[at , a+] = 5ij , [ai ,aj] = [a+ , a+] = 0, (B8) 



and now want to fix the interpretation of the two diagonal 
operators Qo and Q3. We find 

2/3/^ = Q3/h = afai - a^a2 
Q2/h = afai - a3"a3, 
where the number operator 



(B9) 



N ■: 



(BIO) 



and 



2h 

From (B9) we get 



2* 



= ^i9i-ei)+~ipi-pl). (Bii) 



F =-(2^2-2/3). 



(B12) 



Provided we can interpret Q2 as a charge operator this 
is the well-known eight-fold way relation of Gell-Mann 
and Ne'eman between charge, hypercharge and isospin 
[51, 52]. Inserting (B9) in (B6) and rearranging we get 



= -r 



-K^-3fi2-3y2_4/2. (B13) 



The spectrum of the three-dimensional euclidean 
isotropic harmonic oscillator Hamiltonian (B7) is well- 
known [53, 54]. If we assume the standard interpreta- 
tions in (B12) with Q2 as a charge operator we have 
a relation among quantum numbers (^quantum number ~ 
35^perator/^) from which to determine the spectrum of 
M^/fi^, namely 

n = ni+n2+n3, 71^=0,1,2,3,... (B14) 

Now M is hermittean and therefore must be non- 
negative. With K = i, y = 1, ^3 = ±i as for the 
nucleon, the lowest possible value for n is 1 (where 
= 13/4). For the combination y = 1, = ±i, 
with n = 1 we have K{K + 1) + M"^ = 4 and thus 
only two half integer values are possible for K, namely 
K=l, M2 = ^ and = f , = 1, whereas K=^ 
would demand an impossible negative M^. If we keep 
y = 1 and ask for Z3 = ±| the minimum value of n is 2 
and again we get K{K -|- 1) -|- = 4 allowing for the 
two above mentioned values of K. Note that n = 3 does 
not allow 13 = ±1 if we keep y — 1- 

The relation in (B13) can be cast into an Okubo-form 
by choosing a different set of mutually commuting opera- 
tors. We want to replace the three- component of isospin 
by isospin itself. This is possible because 

I' = If + II + li = \iKl + Mi) + II (B15) 
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and [iff + M^,P] = 0. Thus by rearranging (B15) and 
(B13) we get 



3 V A 



■ Sti^ + [Kl + M|) - Sy^ - 4/2 . 

(B16) 



From (B16) we get the following relation among quantum 

numbers 

K{K+l)+M^ = '^{n+lf-:i+{kl+ml)-\y-'-Ai{i+l). 

(B17) 

Here (fc| + ml) should be considered a single quantum 
number. For a given vahie of + jTig) we may group 
the spectrum in (B17) according io n + y = constant and 
get the Okubo structure (26) for the nominator in the 
centrifugal potential as described in the main text. 



Appendix C: Rayleigh-Ritz solution 



and 



Fir, 



-TT J — TT J —TT 



fifmdeide2d9i. 



(C6) 



When the set of expansion functions is orthogonal, 
(C4) imphes 



Ga = 2EFa, 



(C7) 



from which we get a spectrum of A'^ eigenvalues deter- 
mined as the set of components of a vector E generated 
from the eigenvalues of the matrix F~^G, i.e. 



E = ^eig(F-iG). 



(C8) 



We want to find the eigenvalues E of the following 
equation 



-Ae + v]R{e^,e2, 63) = 2ER{eue2, ds), 



(Cl) 



which is the full eq. (22) in the main text. 

In the Rayleigh-Ritz method [55] one expands the 
eigenfunction on an orthogonal set of base functions with 
a set of expansion coefficients, multiply the equation by 
this expansion, integrates over the entire variable volume 
and end up with a matrix problem in the expansion coef- 
ficients from which a set of eigenvalues can be got. Thus 
with the approximation 



N 



(C2) 



Rn = ^ aifi 
1=1 

we have the integral equation 

/ Rn- (-Ae + V)RNdeid92de3 

-TT •/ — TT •/ —TV 

I Rn ■ 2ERNdeidO2d03. (C3) 

J — v 

The counting variable I in (C2) is a suitable ordering of 
the set of tripples p, q, r in (44) or /, m, n in (33) such 
that we expand on an orthogonal set. The eq. (C3) can 
be interpreted as a vector eigenvalue problem, where a 
is a vector, whose elements are the expansion coefficients 
ai. Thus (C3) is equivalent to the eigenvalue problem 

a^Ga = 2Ea'^Fa, (C4) 

where the matrix c;lements of G and F are given by 

Gim = r r r n ■ (-a^ + V)fmde^d92de3 (C5) 

J — TT J —TV J — TT 



The lowest lying eigenvalues will be better and better 
determined for increasing values of A^ in (C2). For the 
base (44) the integrals (C5) and (C6) can be solved ana- 
lytically, and as (44) is an educated guess based on the 
solution of the separable problem (31) it improves the 
convergence in N for the general problem in (Cl) and 
(C7). 

The exact expressions to be used in constructing G and 
F are given below for the base (44). For r > p, u > s and 
q, t> 1 we have the following orthogonality relations 



^ I'pqr I fstu — 

f pgr [61, 62,03) ■ fstu{0l,02, 6s)d6id92d63 

(C9) 



TT TT /• TT 



TT — TT — TT 



Q'K^5ps5qt5ru- 



Here for a convenient notation we have generalized the 
Kronecker delta 



5i3 



The Laplacian yields 



for i = j A i ^ 

for i = 3 A i = 
for i ^ j 



(CIO) 



f 1^ 

< Jpqr I 



1 del 



del 



I fstu > 



{-p - q -r ) ■ Gtt 6ps5gtSru- 



(Cll) 



The matrix elements for the geodetic potential couples 

the individual base functions and follows from a more 
lengthy calculation below yielding the following expres- 
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sion 



< fpqr I ^1 + ^2 + ^3 I fstu > = 

for p. q, r, s,t,u > : 



+(1 - Sp,)6qtSru ■ 67r3 . 4v^^i^ ■ (-1)'^+-' 



— )- 
1 2<;t 

+(1 - • 67r3 • 4(;5^ 

-(1 - 6pu)5qt5rs ■ ■ 4 (/+„f)2 

-(1 - ^rs)5pu^qt • 67r3 • 4^^^^^ 



(-1)9+* 



forp = 0A,S7^0Au>s: 
247r3 [,5,t,5™4,(-l)- - 5qt5rs^{-ir] 

for j3 = A s = : 
^gf^™ • 67r3(27r2 - ^ + ^) 
+(1 - V)^™ • 487r3 • (-1)9+* 

+ (1 - 5ru)5qt ■ ^^TT^'r^^ ■ (-1)'-+". 



(C12) 



Finally the integrals needed for the matrix elements 
of the centrifugal potential can be solved by a change of 
variables. Exploiting the periodicity of the trigonometric 
functions the domain of integration can be selected to 
suit the new set of variables, see fig. 8 and the section 
below on elementary integrals. 



The result is 



1 1 ^ I f 

sin^ 1(^1-^2) ^ sin^ 1(^2-^3) sin^ 1(^3-^1) ' ^ 

= ^'n^[5ps{5r-q,u-tnn{r + q,u + t)- 5r-q,u+tnn{r + q,u- t)) 
+ Sps{- tnn{r -q,u + t) + 6r+q,u+tnn{r -q,u- t)) 

+ Spu{5r-q,s+tnn{r + q,.s -t) - 5r-q,s-tnn{r + q,s + t)) 
+ Spu{-Sr+q,s+tnn{r - q,s-t) + Sr+q,s-tnn{r - q,s + t)) 
+ ^qt{SSp+r,s+unn{r ~p,u- s) + 66p+r,u-snn{r -p,u + s)) 
+ Sqt{S6r-p,s+u'nn{p + r,u- s) + 55r-p,u-snn{jp + r,u + s)) 
+ 5rs{5p+q,u-tnn{p -q,u + t)- (5p+,,„+(nn(p -q,u- t)) 

tnn{p + q,u + t) + 5p^q,u+tnn{p + q,u~ t)) 
+ Sru{Sp+q,s+tnn{p - q,s-t) - 5p+q^s-tnn{p - q,s + t)) 

+ Sru{-Sp-q,s+tnn{p + q,s -t) + 6p-q^s-tnn{p + q,s + t))], (C13) 

I 



where two more shorthand notations have been intro- and 
duced 



for i = j A I 7^ j I * + i I ~ I * ~ i I for i + j = mod 2 

6Sij = {-1 ioT i = -j Ai^O (C14) = |o otherwise. 

otherwise (C15) 
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-t 



FIG. 8: A change of variables from the horizontal/vertical 

(a;, y) to a 45 degrees incUned system of coordinates (u, t) = 

( — - — , — - — ) needed in order to find the matrix elements of 

the centrifugal potential. The seemingly singular denomina- 
tor in the centrifugal potential is then found to be intcgrable, 
see (C16) to (C18). The domain of integration is expanded to 
suit the new set of variables. This is possible because of the 
periodicity of the trigonometric functions such that functional 
values on the hatched area outside the original domain of in- 
tegration [— TT, n] X [— TT, n] arc identical by parallel transport 
from the hatched area within that same area. 



The factor nn originates from the following rule [56] 
sin mx ■ sin nx 



sin^ X 



-dx 



[ {\ m + n \ — \ m — n \)tt for m — n = mod 2 
1 for m — n = 1 mod 2. 

(C16) 

The integrals in (C16) pop up after the aforementioned 
change of variables which exploits the following trigono- 
metric relations 

cos pa; cos ry — cos rx cos py = 
sin nu sin mt + sin mu sin nt, 



u 
and 



x + y 



x-y 



r+p, m = r — p, (C17) 



cos pa; sm qy — sin qx cos py = 

cos nu sin mt — cos mu sin nt, n = p + q, m=p—q. 

(C18) 



1. Elementary integrals for matrix elements in the 
Rayleigh-Ritz method 

We solve here exemplar integrals for the trigonometric 
basis needed to prove the orthogonality relation (C9), 



the expectation value of the geodetic potential (C12) and 
the centrifugal potential (C13). First the orthogonality 
relation for p > 0, g > 0, r > p (p = is left for the 
reader) . With a slight change in notation for the angular 
variables we seek the scalar product 



< f,9 >= 



fgdxidxjdxk 



(C19) 



between base functions (44) like 



fpqr {xi ,Xj,Xk) = eijk cos pxi sin qxj cos rxk and 

9stu{xuXm,Xn) = Clmn COS pXl sin qXm COS TXn- (C20) 



The three-dimensional integral in (C19) factorizes into 
three one-dimensional integrals and using the orthogonal- 
ity of the trigonometric functions on [— 7r,7r] we readily 
have 



^ fpqr-; 9stu ^ — (^j>5 

which for p > 0, q > 0, r > p and s>0,t>0,u>s 
reduces to (C9). To obtain the expectation value (C12) 
of the geodetic potential we use the same kind of fac- 
torization together with the following list of elementary 
integrals 



I X^ cos px cos Sxds = {—lY'^^^-K-—^ p s 

J-T, (" 



(p2 — s^) 



J — TT 



47r 



cossa;rfa; = ^(— 1)*, p = 0,s^O 



2 2,^ TT 

X COS pxax = 1 K , p = s 

3 2p2' ^ 



/7r 
-TT 



2 . 2 

X sm qxax = — — — ^ , q = t 
3 zq-' 



/ x^singa; sin ta;rfa; =(—1 )*"'"* 47r 

J — TT 

/TT 
x^ COS pa; sin fa; dx = 0. 
-TT 



2qt 



(g2 - t^y 



(C22) 
(C23) 

(C24) 

(C25) 

\l\ ^ \t\ 
(C26) 

(C27) 



For the Laplacian and for the centrifugal potential we 

make another slight change in notation for our angular 
variables and rewrite our base functions (44) as a sum of 
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subdeterminants 



/pgr (2^! 2/) 



COS pa; cos py cos pz 
sin ga; sin qy sin 
cosra; cosry cos rz 



= cos pz 



— sm qz 



- cos rz 



sm qx sin g'?/ 

cosra; cosry 

cos pa; cos py 

cosra; cosry 

cos pa; cos py 

sin ga; sin 



Prom (C28) we then get terms like 



fpqr{x,y,z) = -p^COSpZ 



(C28) 



sm qx sm 
cos ra; cos ry 



+q sinqz 



-r cos rz 



cos pa; cos py 
sin qx sin 

cos pa; cos py 
sin ga; sin qy 



(C29) 



for the Laplacian. Again the integral for the expectation 
value factorizcs into one-dimensional integrals where the 
orthogonality of the trigonometric functions can be ex- 
ploited to obtain (Cll). 



For the centrifugal term the three-dimensional integral 
does not readily factorize. We need a change of variables 
which suits the mixing of the variables in the denomi- 
nators. Due to the arbitrary labelling of our angles we 
have 



< /p. 



sin 



2 1 



{x-y) 



+ 



1 



sm 



2 1 



3 < fpqr \ ^ . 2 1 



1 



+ 



sin 



2 1 



{z-x) 



I fstu > 



sin Ux-y) 



I fstu > ■ 



With 



fpqr fstu 



COS pa; COS py cos pz 
sin qx sin qy sin qz 
cosrx cosry cos rz 



cos sx COS sy cos sz 
sin tx sin ty sin tz 
cosux cos uy cos uz 



(C30) 



(C31) 



we can use the subdcterminant expressions in (C28) to get e.g. a factor SSqiir from the z-integration of the term 
involving the two sines while the product of the two corresponding subdeterminants is used for a shift of variables, 
see below. We have 

cos pa; cos py cossx cos sy 
cosrx cosry cosux cos uy 
= (cos pa; cos ry — cos rx cos py) (cos sx cos uy — cos ux cos sy) 
_ 11, 
22 



[cos(pa; — ry) + cos{px + ry) — cos{rx — py) — cos{rx + py)] 



•[(cos(sa; — uy) + cos{sx + uy) — cos{ux — sy) — cos{ux — sy)] 
= [— sin ^ (pa; — ry + rx — py) sin ^ (pa; — ry — {rx — py)) 



-sin - (pa; 



ry + rx +py) sin ^(pa; + ry - {rx +py))] 



■ [— sin ^ {sx — uy + ux — sy) sin ^ {sx — uy — {ux — sy)) 

— sin -(sa; + uy + ux + sy) sin -^{sx + uy — {ux + sy))] 

= [sin nv sin mw — sin nw sin mv] [sin kv sin Iw — sin kw sin Iv] , 



where 



p + r, m = p ■ 



r, k = s + u, I = s — u, w = — - — and v = — - — . 



(C32) 

(C33) 
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Since both nominator and denominator in (C30) are 
trigonometric functions we can exploit their periodicity 
to enlarge the domain of integration and make a shift of 
variables to w and v. An integration over the original 
domain is namely half the value of an integration over 
the enlarged domain in fig. 8, thus 



/ dxdy = ^ dw' / dv' 

/2'K p2'K /"TT /'TT 

dw' / dv' = dw dv, 
-2-K J — 2-K J — TT J —TT 



1 1 
22 



where 



y 



and 



x + y. 



(C34) 



The factor ^ in the second expression is just from the 
change of coordinates dw' du' =^ 2dxdy and the domain 
of integration is still not enlarged, but limited by piece- 
wise linear functions v[ and V2 ■ In the third expression 
then we double the area of integration to lift the coupling 
between w' and v'. In the last expression we just rescale 
our variables to suit our needs in (C32). 

With the coordinate transformations in (C34) we can 
use (CI 6) to get the final result 



3 < /pgr 1.21 



sin l{x-y) 



fstu ^ 



„ _ , , , [sin nv sin mw — sin nw sin mv] [sin kv sin Iw — sin kw sin Iv] 
3ogt • TT / dw / dv 



TT J — TT 



sm V 



= 3Sqt ■ TT[5m,l • 7r(| n + fc | - | n - fc |) • TT - 5mk ■ nil n + 1 \ - \ n - I \) ■ TT 



■ Sni ■ tt{\ m + k \ - \ m - k \) ■ TT + Snk ■ tt{\ m + I \ - \ m - I \) ■ tt], 



(C35) 



which is a specific example of the general result (C13). we will project on ^ (1 — x)T in order to probe on 

xPij,. With a probability amplitude interpretation of $ 
we then have 

fT{x)-dx = {dRu=c^p(eiT)f-dO, where 6 = tt£_, (Di) 

fxix) is the sought for distribution function and R{0) 
is the proton analogue of (44) with the period doublings 
in (41) mimicked by making q half integer, see (Dll) 
below. The distributions in fig. 7 are traced out with the 
generators 



Appendix D: Toroidal parton distributions 

Inspired by Bettini's elegant treatment of parton scat- 
tering [57] we generate distribution functions via our 
exterior derivative (9). In short the derivation runs 
like this (with H = c = 1): Imagine a proton at rest 
with four-momentum P = {0,Eo). We boost it to en- 
ergy E by impacting upon it a massless four-momentum 
q — (q, -E — Eq) which we assume to hit a parton xP. 
After impact the parton carries a mass xE. Thus 



T 



(xP^ + q^)-ixP'' + q'')=x^E^, 
from which we get the parton momentum fraction 

2Eo 

E + Eq 

Now introduce a boost parameter 

^ E-Eo 2-2x 



E 



2-x 



(Dl) 



(D2) 



(D3) 




and Td 




With (D3) in (D4) we have in general 

X ■ fT(x) = X ■ dP? ,2-2x . 



TT-2 



'^T) {2-xf 



(D5) 



(D6) 



which we shall use to track out an orbit on the u(3) 
torus (^ = 1 corresponds to a; = and vice versa). 
With the toroidal generator T as introtangling momen- 
tum operator we namely have the qualitative correspon- 
dence q^ E - Eq [l- x)E ~ (1 - x)T. That is. 



The peak position of the distribution functions gene- 
rated by (D6) depends on the 'direction' of T relative to 
(Ti,T2,ri). The choices (D5) are 'skew' fractionations 
of the charge operator Q2 in (Bll) which in a matrix 
representation reads 



Q2/n = 




(D7) 
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Since the tracking in (D6) is done from laboratory 
space we project on a fixed basis like in (A8) when we 
calculate the exterior derivative dR. Any toroidal track- 
ing generator T is a tangent vector and can be expressed 
in a fixed basis of the tangent space of the torus (the 
Cartan subalgebra). We write 



T = aiTi + a2T2 + asTs, 



(D8) 



where aj are the coefficients of expansion and the T^'s 
are the toroidal generators also used in (4). In quark 
language in the end we want to project dR on all three 
colour degrees of freedom - cf. (A4) - and sum the contri- 
butions to got the total amplitude for the specific flavour 
amplitude along the specific tracing generator. In general 
we have [14, 58, 59, 60, 61] (see appendix E) 

3 3 

J2dRu{dj) = Y,dRu{mTj) 

d 



J2dRu{de^p\^,^,^,^^ig^)) 



(D9) 



where u = exp(6li?Ti + O^iT^ + 9:^iT:i) and R* is the pull- 
back of R to parameter space from the manifold. We 
shall immediately omit the asterix on R* again and we 
have finally 



{di + d2 + ds) 

dR dR dR' 



dO-i 802 863 J Uei.«2,e3) = (e-ai,e-a2,0-a3) 



D{e-ai,e-a2,e-a3) = DT{e), 



(DIO) 



where wc have introduced a shorthand notation for the 
directional derivative. 

In our specific case in fig. 7 we use a first order term 
(C2) from (44) with {p, q, r) = (0, i, 1) 



Vi,i(^i'^2,e3) = 



N 



1 



1 



1 



sin 2 
cos^i 



\9i sin 5 6*2 sin i 
= ^- cos 62 cos 6z 



, (Dll) 



where A'^ is a normalization constant. For the particular 
case of (Dll) we get a specific expression for the function 
D introduced in (DIO) 



ND{ei,e2M = 

— - COS — • (cos 6/3 — COS 62) — sm 0\ ■ (sm — sm.—) 

, ■'■ ^2 / a a \ , ■ a I ■ • ^1 ^ 

2 ~2 ' ^ ~ cosWi) -|- SU1O2 ■ (sm — sm — ) 

— - cos — • (cos O2 — COS Oi) — sm O3 ■ (sm — sm — ). 

(D12) 

And tracking with r„ from (D5) we get the final result 
for our distribution function 



Xfu{x) 



D{e --,9 -0,6 ■{-!)) 



tdet 

' dx ' 



where 6 = n^. 



(D13) 



Expressed in the parton fraction x (D13) reads 



^ , 2 - 2a; 2 2 - 2x' „ 2 - 

^ ■ TT)""^; • t),7r— - 

2-x 3 2-x 2 



2x 



(-1) 



TT- 2 



(3 



(D14) 



r 



The result (D14) is shown in fig. 7 together with the 
result of tracking with Trf. Note that a:- /tm (a;) |macc~ 0.63 
and X ■ ,fTd{x) \max'^ 0-32 in rather fine agreement with 
u„ and dt, in the distribution functions shown in fig. 7, 
and which arc cxpcirimcntally extracted via NLO-QCD 
parametrizations. Note also that 



' xfTu{x)dx / / xfTd{x)dx 
Jo 

0.2722/0.1437 = 1.89 « 2, 



(D15) 



which for a first order approximation (Dll) is in rather 



close agreement with the relative constituent quark 
flavour niodcil content of the proton according to the stan- 
dard model, namely 2u:ld. Actually the ratio between 
the distributions themselves yields 



fTu{x)dx / I /Td(a;)da; = 1.33.../0.66 = 1.996.... 



(D16) 

We should also note that when comparing with estab- 
lished distribution functions one should be aware that 
these are slightly scale dependent as seen when compar- 
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ing those for 10 GcV^ with those for 10,000 GeV^ from 
the particle data group [62]. One should expect closer 
correspondence the lower the energy scale. Actually we 
would like to see data extracted at « GeV^, which 
would represent the purely allospatially derived distribu- 
tion functions. Luckily though, it seems that 10 GeV^ is 
small enough to make the bridge. We could call it the 
bridge of soft 'deep' inelastic scattering. The softer the 
momentum transfer gets the more terms one needs in the 
QCD-expansions through which the data are currently 
analyzed. The distributions in the insert in fig. 7 are 
paramctrizations from next to leading order, NLO calcu- 
lations. What we would need in order to reach allospa- 
tial dynamics would be N°°LO-QCD at = GeV^. 
This would be the confinement domain which we have 
attacked in the present work 'from within', but which is 
the domain where pertubative QCD is expected to break 
down, i.e. where N"LO-QCD is not expected to converge 
for n — )• 00. 



Appendix E: Vector fields, derivations and forms on 
smooth manifolds 

To understand the use in (A6) of left invariance of the 
coordinate fields and linearity of the coordinate forms 
actually requires a careful application of the different ex- 
pressions for the directional derivative. And to arrive 
from the first expression in (D9) to the last expression of 
the same equation specifically requires a consistent def- 
inition of differentiation on smooth manifolds. In the 
present appendix we shall outline the basic concepts. The 
appendix is based almost exclusively on [58]. We have 
only added two corollaries of particular relevance in the 
present context. Here it comes: 

Let Z be the generator of a directional derivative of ^ 
at a certain point u in our manifold u(3). Repeating eq. 8 
of [14] we then have the following equivalent expressions 
for the directional derivative 

[*] = Z [*](u) = (d*)„(Z) = ^\I/(«exp(fZ)) |t=o. 

(El) 

The identity between the second and the third expression 
in (El) is just a question of notation. 

Referring to fig. 9 we namely have for the directional 
derivative of real functions / on a smooth manifold M 
with tangent space TM the following 

Definition 1: 

Let / : M i? be a smooth function and Xp £ TMp. 

The number dfp{Xp) is called the (directional) derivative 

of /in the direction Xp and is also denoted Xp{f). 

(E2) 

The directional derivative Xp{f) fulfils the following 



identities 

Xpif + g) = Xp{f) + Xpig) f,g&C^{M) (E3) 

X(A • /) = A • Xp{f) A € i? (E4) 

Xpif ■ g) = Xpif) ■ g{p) + f{p) ■ Xp{g). (E5) 

where f + g and / • g are the functions given by addi- 
tion and multiplication of functional values. We see that 
the directional derivative applies and acts algebraically 
in analogy to the usual act of differentiation. 
We now state the following 

Definition 2: 
A map Dp : C°°(M) ^ R which fulfils (E3)-(E5) 

is called a derivation at the point p. (E6) 

We see that any tangent vector gives rise to a derivation. 
It is remarkable and important for our purpose, that the 
opposite is also true: 

Theorem, 1: 
Any derivation Dp at the point p 

is a directional derivative. (E7) 

This theorem is used in order to write the identity be- 
tween the third and the fourth (last) expression in (El). 
We shall leave out the proof of (E7) which needs some 
lemmas that would lead us astray in the present con- 
text. The conclusion is that there exists a 1-1 correspon- 
dence between tangent vectors Xp e TMp and deriva- 
tions Dp : C°°{M) R aX the point p: The derivation 
corresponding to Xp is 

Xp{f) = dfp{Xp) or (E8) 
Xp{f) = Y,a,^\p, (E9) 

i * 

where Xp = I], a^^jp and §fjp= ^^^^\x(p) for an 
atlas {U,x), that covers the manifold M. 

The idea is illustrated in fig. 9 locally for a partic- 
ular local map in the book of maps, the atlas, which 
consist of enough maps to parametrize all of M via lo- 
cal maps that are joined smoothly in their overlapping 
regions. The manifold M is depictured in fig. 9 as a 
surface in an embedding space of specific dimension. So 
now in short we have a global definition of derivation and 
can faithfully derive locally by pulling back the derivation 
on the manifold to ordinary derivations in a parameter 
space R™ of the same dimension m as the manifold M. 
Note that the presence of the spikes in the parametric 
potential (17) shown in fig. 2 does not signify that the 
geodetic distance (3) used for our potential (2) in our 
defining Hamiltonian (3) would not be smooth. Quite 
the contrary: The geodetic distance, and thus its square, 
are smooth fimctions on all of u(3). The spikes in fig. 
2 only indicate that one cannot define a single map to 
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FIG. 9: Derivation of a real function / : M — > R at point p 
in the mainfold M is defined by using a local smooth map x : 
M — >■ _R™ to pull back the problem to an ordinary derivation 
on i?™ by using the puUback function / o : ET^ — >■ 7?. 
One can differentiate f o in the ordinary way. This idea 
is readily generalized to a complex valued function and in our 
case the manifold M could be u(3) and the complex valued 
function / could be either the wavefunction or its measure 
scaled partner 3>. 



cover u(3) globally ("You cannot peel an orange without 
breaking the skin"). This is also why we have to limit 
the expression (2) to the interval — tt < Oj < tt for our 
parametrization which we need to make our Hamiltonian 
operational for actual calculations in (16). But luckily 
we can use a common parametrization of the eigenvalues 
of u € w(3). 

To proceed further we need the following 

Statement 1: 
The smooth vector fields T{M) on M constitute 
an (infinite dimensional) real vector space. (ElO) 

My personal notes in the lecture notes [63] explains the 
infinity with these comments: "There is an infinity of 
atlases and infinitely many coefficients to At any 

point p G M there is an m-dimensional infinity of tan- 
gent vectors, (and) these can be combined in infinitely 
many ways at diff'erent p's." To understand the statement 
(ElO) wc need some more math. First let M™ C R'' be 
a smooth manifold. We shall then start out with the 
following [64] 

Definition 3: 

A function X : B!" which fulfils X{p) e TMp 

for all p e M is called a vector field on M. If M is 
continuous (smooth), the vector field is said to be 

continuous (smooth). (Ell) 

Please allow me to cite my teacher's comment on the 
above definition: "Properties of vector fields and their 
dynamical systems on manifolds are studied locally by 
translating the previously treated euclidcan case by the 
help of atlases." This is a key point in our application in 
(D9) of the beautiful mathematical machinery exposed 
in this appendix. Let us now continue with his explana- 
tion following my citation: Let ([/, x) be a chart on M™. 



Corresponding to this we have an induced base for any 
point p G U: 

^\p=d{x-%^p){ei), (E12) 

where d{x-^)^^p) : i?™ TM™ C R'' and where e, is 
the i'th standard unit vector in R^. Since ei, . . . , is 
a base for i?™ 

^^^^^ 

becomes a basis for TM'^ . I made a remark in his notes 
that "d(x~^) is injective since a; is a diffeomorphism" (per 
definition of the smoothness of M, OT 2012). Perhaps 
we should have mentioned earlier, that U C M is an open 
neighbourhood of p G U, which is where we have defined 
(directional) derivatives above (E2). 

Going back through the necessary mathematical con- 
cepts for our purpose, we give for the sake of complete- 
ness: 

Definition 4 '■ 

Let p e M" C i?*" be a point on a smooth m-dimensional 
manifold and {U, x) a chart around p with x{p) = a. 
The tangent space TMp for M"* in the point p is the 
image of the linear map 

rf(a;-i)„ : i?™ i?^ (E14) 



Definition 5: 

Let Af" G R^ and A^" € i?' be smooth manifolds and 
/ : M™ ^ iV" a smooth function. By the differential of 
/ in the point p is understood the linear map 

dfp : TMp ^ TNf(p) (E15) 
defined in the following way: Let F be a neighbourhood 
of p in R'^ and / : F ^ i?' a smooth continuation of / 
(/ I y n M = /) (or in words: / 's restriction to M 
is /, OT 2012). Then for ^ e TMp the diffential is 

dm) = djpii). (E16) 

We are now ready to explain statement 1 in (ElO), 
which we shall need for essential algebraic manipulations 
below. Statement 1 might be obvious in connection with 
the lecture notes, but here we need to strengthen the 
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statement by formulating the following 
Corollary 1: 

Let df : TM"" ~>TR = Rhe the diffential (the exterior 
derivative) of a function/ : M™ — )• R. Further let 

^j-(/) = E«^^lf (E17) 

be (directional) derivatives of / in the directions Xp 
respectively Yp. Then 

^pif) = E '^'a^lf ' where Cj = a, + 6^, (E19) 

i * 

will be a directional derivative in the direction Zp. 
In other words 
{Xp + Yp){f) = Xpif) + Yp{f). (E20) 

Proof of corollary 1: 
It should be clear that TM^ is a linear vector space and 
thus it is clear that Zp e TMp. We also remarked in the 
text following (E6), tliat Zp will be a derivation. This is 
obvioiis from definition 1 in (E2) when comparing with 
definition 2 in (E6). Now from theorem 1 in (E7) we 
then have that Zp will also be a (directional) derivative. 
Q.E.D. It shall be essential to us that Zp in (E19) is ex- 
panded on the same basis as that of Xp and Yp . This is 
possible because one can use the same smooth continu- 
ation / for all the three derivatives Xp, Yp and Zp. We 
need a second 

Corollary 2: 

Let / ■.M'^-^R then 
df{X + Y)=df{X) + df{Y), (E21) 
where X, y e TM" 
and df{XX) = A • df{X), (E22) 
where X G R. 

Here we have suppressed the index p, but (E21) and 

(E22) should be understood pointwisc. It should be ob- 
vious that (E21) in corollary 2 is just a formal way of 
expressing the statement in definition 5, that dfp is a lin- 
ear map, and that (E22) in corollary 2 is just another 
way of stating one of the characteristics of derivatives, 
namely (E4). However, we have not proven the linearity 
of forms, which has been the formulation (exterior deriva- 
tive) primarily used in the main text of the present work. 
Thus we wanted to state this small 

Proof of corollary 2: 

df,iX + Y) = iXp + Yj,){f) = 
Xpif) + Ypif) = dfp{X) + dfp{Y). (E23) 



(E2), in the third expression wc used (E20) and in the 
last expression we returned to the language of forms via 
(E2). To prove (E22) we write 

d/p(AX) = XXp{f) = Xdfp{X). (E24) 

In the second expression wc used again definition 1 in 
(E2). The last expression finishes the proof. Q.E.D 

Of course the "Proof of corollary 2" is not really a 
proof, there was nothing new to prove. The 'proof is only 
meant to show how the algebra works when one switches 
between the different notations Xp{f) and dfp{X) for the 
directional derivative. 

Corollaries 1 and 2 arc readily generalized to / : 
M"* — C and A G C. With these two corollaries we 
now have our mathematical machinery in shape to pro- 
ceed from the first to the last expression in (DIO) in order 
to calculate (D12). 

The linearity of forms expressed in corollary 2 by the 
eqs. (E21) and (E22) together with the generalization 
(E2) of derivation in definition 1 and finally the existence 
of induced bases like in (E12) pull back the problem to 
i?™ so that the partial derivatives pull back to give results 
like (D12). We have thus presented the derivation on 
the smooth manifold u(3) in an operational form which 
makes possible actual calculations. 



In the second expression we used the equivalent notation 
for directional derivatives introduced in definition 1 in 
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Tables 



TABLE I: Comparison of the 10 lowest lying eigenvalues E of (l)/(22) for isospin | and hypercharge 1. The eigenvalues are 

given without and with the curvature and centrifugal potentials from a Raylcigh-Ritz calculation with 1800 base functions of 
the type (44) corresponding to p = 0, 1, 2, 14; g = 1, 2, 15; r = p + 1, 15 . Minimum values have been used for K"^ and 
M'^, namely K{K + 1) = |(| + 1) and = ^ , following the interpretation in (B14). The corresponding shifts in energy 
between calculations without and with the curvature and centrifugal potentials lie within 50 MeV as claimed in the main text. 
The somewhat curious order in the list of the alleged candidates follows from choosing them according to the eigenvalue of the 
charged partners from approximate solutions in sect. 4. The value 1723 MeV in brackets in the second column is the eigenvalue 
of such a charged partner. 



Eigenvalues 


E 


Eigenvalues 


E 


Candidate 


Relative 


no curvature 




curvature 






shift from 


no centrifugal (31) 


MeV 


centrifugal (22) 


MeV 




(31) to (22) 


4.47 


939.6 (fit) 


4.38 


939.6 (fit) 


n, p 


-2.0 


6.22 


1308 


6.10 


1309 


N(1440) 


-1.9 


6.57 


1381 


6.53 


1401 


N(1535) 


-0.6 


8.19 


1722 


8.12 


1742 


N(1675) 


-0.9 


8.32 


1749 


8.42 


1806 


N(1720) 


1.2 


8.39 


1764 


8.53 


1830 


N(1680) 


1.7 


9.17 


1928 (1723) 


9.35 


2006 


N(1650) 


2.0 


10.29 


2163 


10.37 


2225 


N(2250) 


0.8 


10.36 


2178 


10.54 


2261 


N(2190) 


1.8 


10.48 


2203 


10.60 


2274 


N(2220) 


1.1 
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TABLE II: Eigenvalues of the parametric group space chopped harmonic oscillator Schrodinger equation (32) calculated with 
1500 collocation points. The same results are obtained by 1) iterative integration, 2) McLaurin series [65] and 3) collocation, 
see table 6. The mutual discrepancies are at the level of 10~*. This fine agreement among the different methods lends support 
to the Rayleigh-Ritz method also for solving the full eq. (22). Note that the lowest eigenvalues as expected are close to 
those of the ordinary harmonic oscillator since the lowest states live in the neighbourhood of origo and thus do not 'feel' that 
the potential is chopped into periodicity. Moving up to higher levels the eigenvalues differ more and more from those of the 
harmonic oscillator as indicated in fig. 5. 



Level Eigenvalue Diminished Augmented 

1 0.4998047079793375 0.5001727903915900 

2 1.502988968183189 1.496433950157817 

3 2.471378779213570 2.522629649224744 

4 3.600509000413400 3.377236031951678 

5 4.218515963091988 4.803947526779894 

6 6.197629004032325 5.160535373287425 

7 6.383117406428158 7.820486992163699 

8 9.688466291114409 7.922699153838795 

9 9.751335596178837 11.80644675634627 

10 14.17552754458349 11.84897047349423 

11 14.20637080732072 16.79575229446833 



TABLE III: Scarce singlet states. Eigenvalues based on Slater determinants of three cosines up to order 20 analogous to (44). 
The first column shows eigenvalues of the approximate eq. (31) and the third column shows eigenvalues of the exact eq. (22). A 
singlet 579-likc resonance is predicted at 4499 MeV in the free charm system (2455)1)^ slightly above its threshold at 4324 
MeV. It should be visible in neutron diffraction dissociation experiments like those in ref. [37] . The rest masses are predicted 
with the common fit of table 1 where 939.6 MeV corresponds to the ground state N. 



Singlet Toroidal Singlet Toroidal Rest mass 
approximate (31) label exact (22) partner MeV/c^ 



7.18957249 


1 3 5 


7.12174265 


1 3 5 


1526 


9.35680923 


1 3 7 


9.57104964 


13 7 


2051 


11.11924735 


1 5 7 


11.29403818 


1 5 7 


2420 


12.71754526 


1 3 9 


13.25048053 


1 3 9 


2839 


13.09274266 


3 5 7 


13.28113847 




2846 


14.44940740 


1 5 9 


14.96408154 


15 9 


3206 


16.40861790 


3 5 9 


16.92132768 




3626 


16.66054860 


1 7 9 


17.30060283 


1 7 9 


3707 


17.17694213 


1 3 11 


18.00904018 


1 3 11 


3859 


18.63196975 


3 7 9 


19.25767851 




4126 


18.92139903 


1 5 11 


19.73271915 


1 5 11 


4228 


20.37744401 


5 7 9 


20.99400214 


5 79 


4499 


20.89101474 


3 5 11 


21.71097283 


3 5 11 


4652 


21.07660092 


1 7 11 


22.04093052 


1 7 11 


4723 
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TABLE IV: Consistency of results from expanding on Slater determinants of 1-dimensional functions like (33) and (44) to solve 
the full cq. (22). In the first column are the lowest lying N-states constructed from two cosines and one sine as in (44). In 
the second and third columns arc the lowest lying states constructed from parametric oigonfunctions as in (33). Note that the 
states based on parametric eigenfunctions cover the N-statos fully and further include also states that are identified as either 
singlet or A-states. The parametric base eigenvalues use the eleven lowest lying solutions to (32) with respectively 47^ and 29^ 
base points for each function. The rest masses in the last column are predicted from a fit to the ground state of the parametric 
base where 939.6 MeV corresponds to the ground state eigenvalue 4.364. 



Trigonometric Parametric base Toroidal Rest mass 
base eigenvalues labels 
N-states 47^ 29^ Mev/c^ 



4.385 


4.364 


4.345 




939.6 (fit) 




5.568 


5.546 A-state 


124 


1199 


6.103 


6.087 


6.067 






6.537 


6.517 


6.493 








7.054 


7.033 Singlet 


135 


1519 




7.538 


7.508 A-state 


234 


1623 


8.119 


8.087 


8.060 






8.424 


8.368 


8.333 






8.533 


8.513 


8.476 






9.356 


9.326 


9.290 
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TABLE V: Spin and parity conjecture on allospatial states compared with the series of all the observed 4-star N-resonances. 
The assignments are only tentative. An alternative assignment around the singlet state 137 is shown in brackets. 



{-!)"''■ ; I, m,n L2i,2j Name RPP rating 



+ 123 


1/2+ 


Pll 


n, p 


+ * + =t: 


+ 125 


1/2+ 


Pll 


N(1440) 


**** 


-134 


1/2- 


Sll 


N(1535) 


**** 


+(135) 


1/2- 


Sll 


N(1650) 


**** 


- 235 


3/2- 


D13 


N(1520) 


**** 


-hl45 


3/2+ 


P13 


N(1720) 


**** 


+ 127 


5/2+ 


F15 


N(1680) 


**** 


-136 


5/2- 


D15 


N(1675) 


**** 


+ (137) 






N(2040) 




+ 345 


(9/2-) (G19) (N(2250)) 




- 237 


7/2- 


G17 


N(2190) 


**** 


+ 147 


9/2+ 


H19 


N(2220) 


**** 


- 156 


9/2- 


G19 


N(2250) 


**** 



-257 



11/2- nil N(2600) 



*** 
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TABLE VI: Comparison of numerical results for the eigenvalue of the ground state. The seperable problem (31) has been 
solved by four different methods three of which gives a set of eigenvalues for the one-dimensional problem (32) from which the 
eigenvalues for the three-dimensional problem (31) is constructed. These eigenvalues can be used to check the Rayleigh-Ritz 

method for solving the three-dimensional problem directly. Mutual discrepancies arc due to the finite expansions in the different 
methods. The fine agreement among the different methods lends support to the Rayleigh-Ritz method also for solving the full 
eq. (22). 



ID-level Iterative MacLaurin Rayleigh-Ritz Collocation Collocation 

number integration series [65] 1800 1500 1500 

base functions points points 

Comal - MathCad Matlab m„ Matlab rUp 



1 0.499804708 0.499804704 - 0.499804708 0.500172790 

2 1.502988981 1.502988968 - 1.502988968 1.496433950 

3 2.471378882 2.471378899 - 2.471378779 2.471378779 
Sum 4.474172571 4.474172571 4.47417271 4.474172455 4.467985519 
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